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Abstract. In this communication, we show how asteroids observations from the Gaia mission can be 
used to perform local tests of General Relativity (GR). This ESA mission, launched in December 2013, will 
observe -in addition to the stars- a large number of small Solar System Objects (SSOs) with unprecedented 
astrometric precision. Indeed, it is expected that about 360,000 asteroids will be observed with a nominal 
sub-mas precision. 

Here, we show how these observations can be used to constrain some extensions to General Relativity. 
We present results of SSOs simulations that take into account the time sequences over 5 years and geometry 
of the observations that are particular to Gaia. We present a sensitivity study on various GR extensions 
and dynamical parameters including: the Sun quadrupolar moment J 2 , the parametrized post-Newtonian 
parameter /3, the Nordtvedt parameter 77 , the fifth force formalism, the Lense-Thirring effect, a temporal 
variation of the gravitational parameter GMq (a linear variation as well as a periodic variation), the Standard 
Model Extension formalism, ... Some implications for planetary ephemerides analysis are also briefly 
discussed. 
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1 Introduction 

This year marks the centenary of the publication of the classical theory of General Relativity (GR), the current 
paradigm to describe the gravitational interaction. Since its publication in 1915, GR has been conhrmed by 
experimental observations. Although very successful so far, it is nowadays commonly admitted that GR is not 
the ultimate theory of gravitation. Attempts to develop a quantum theory of gravitation or to unify gravitation 
with the others fundamental interactions lead to deviation from GR. Moreover, observations requiring the 
introduction of Dark Matter and Dark Energy are sometimes interpreted as a hint that gravitation presents 
some deviations from GR at large scales. 

GR is built upon two fundamental principles. The first principle is the Einstein Equivalence Principle (EEP) 
which gives to gravitation a geometric nature. More precisely, EEP implies that gravitation can be identified 
to space-time curvature which is mathematically described by a space-time metric If the EEP postulates 
the existence of a metric, the second building of GR specifies the form of this metric. The metric tensor is 
determined by solving the Einstein field equations which can be derived from the Einstein-Hilbert action. 

The EEP has been thoroughly tested (Will 2014, and references therein) by considering Universality of Free 
Fall experiments (Schlamniinger ct al. 2008; Adclberger et al. 2009; Williams et al. 2009; Miillcr et al. 2012), 
constancy of the constants of Nature (Uzaii 2011; Rosenband et al. 2008; Guena et al. 2012; Miiiazzoli & Hecs 
2014; Hecs et al. 2014e), redshift experiments (Vcssot et al. 1980; Delva et al. 2015), anisotropy in the velocity 
of light (Whll 2014), ... 

On the other hand, in the last 50 years, mainly two formalisms have been used to test the form of the metric: 
the parametrized post-Newtonian (PPN) formalism fully described in Will (1993) and the fifth force formalism 
described in Talmadge et al. (1988); Adclberger et al. (2009). The PPN formalism is making a nice interface 
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between theoretical developments and observations by parametrizing deviations from GR at the level of the 
space-time metric by means of 10 dimensionless coefficients. The fifth force formalism considers a modification 
of the Newtonian potential of the form of a Yukawa potential. The parameters entering these formalisms 
have been severely constrained by diverse measurements (for a review, see Will 2014): spacecraft tracking 
(Bcrtotti ct al. 2003; Konopliv ct al. 2011), Lunar Laser Ranging (Williams ct al. 2009; Miillcr ct al. 2012), 
planetary ephemerides (Ficnga et al. 2011; Pitjeva & Pitjev 2013; Verma et al. 2014; Ficnga ct al. 2015), Very 
Long Baseline Interferometry (Lambert & Lc Poncin-Lafitte 2009), etc... Despite very stringent constraints, 
we have strong theoretical motivations to pursue these tests and to consider formalisms beyond the standard 
PPN and fifth force formalism (for some motivations, see for example Hees et al. 2012, 2014c, d). 

Launched in December 2013, the ESA Gaia mission is scanning regularly the whole celestial sphere once 
every 6 months providing high precision astrometric data for a huge number (r; 1 billion) of celestial bodies. 
In addition to stars, it is also observing SSOs, in particular asteroids. One can estimate that about 360 000 
asteroids will be regularly observed. The high precision astrometry (at sub-mas level) will allow us to perform 
competitive tests of gravitation and to provide new constraints on alternative theories of gravitation. These 
constraints will be complementary to the ones existing currently since relying on different bodies, on different 
type of observations and therefore sensitive to other systematics. 

In this communication, we report the first results of a sensitivity study of Gaia SSOs observations to several 
modifications of the gravitation theory. Gorrelations between the different dynamical parameters (in particular 
with the Sun quadrupolar moment) are also assessed. 


2 Simulations and parameters estimation 

We have simulated the trajectories of asteroids using the standard post-Newtonian equations of motion in a 
heliocentric frame. The mutual interactions between the asteroids are neglected, the Sun oblateness parameter 
J 2 is considered and the perturbations from the different planets and the Moon are modelled by using INPOPlOe 
ephemerides (Ficnga ct al. 2013). The initial conditions used for the simulations are provided by the ASTORB 
database* and a match between the asteroids trajectories with the Gaia scanning law is performed to find the 
observation times for each asteroid. Simultaneously with the equations of motion, we integrate the variational 
equations (for a detailed presentation of the method, see Hcstroffer et al. (2010); Mourct & Mignard (2011); 
Mouret (2011)). The simulated asteroid trajectories are transformed into astrometric observables as well as 
their partial derivatives with respect to the parameters considered in the covariance analysis. The parameters 
considered here are twofolds: (i) local parameters specific to each asteroid (e.g. their 6 initial conditions) and 
(ii) global parameters that are common for all asteroids (the Sun J 2 and the parameters that characterize 
the gravitation theory). Our sensitivity study is performed by computing the covariance matrix to assess the 
correlations between the estimated parameters and the formal uncertainties reachable using Gaia observations. 
In this proceedings, the simulations performed include only 10,000 asteroids and the astrometric accuracy used 
is 0.2 mas. Full simulations of 360,000 asteroids with updated and better precision are ongoing and will be 
published in the future. 


3 Sensitivity study to various global parameters 

In all the following, we consider modifications of the gravitation theory that do not produce any important effects 
on the light propagation. Therefore, we concentrate on the orbital dynamics only and neglect any modification in 
the propagation of the light. Simulations using the Time Transfer Formalism (Teyssandicr & Lc Poncin-Lafittc 
2008; Hees et al. 2012, 2014b, a) have been performed to ensure that the effects of the considered gravitation 
modifications on the light propagation can safely be neglected. 

In the following subsections, we report a sensitivity study performed by a global inversion that includes the 6 
initial conditions for each of the 10,000 asteroids, the Sun J 2 and the parameters characterizing the gravitation 
theory. 


see http://www.naic.edu/~nolciii/astorb.html 
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3.1 Sun quadrupolar moment, PPN parameter and Nordtvedt effect 

First of all, we consider the standard PPN parameter fd (Will 1993, 2014). It is well known that this parameter 
is highly correlated with the Sun quadrupole moment J 2 . Indeed, the secular advances of the perihelia produced 
by these two parameters are given by 



(2 + 27 - /3)n 


GM 

0^0(1 — e^) 


3 J 2 RI 
2^a^{l-e^Y ’ 


(3.1) 


where lo is the argument of the perihelia, GM is the Sun gravitational parameter, c the speed of light, 7 and 
j3 are the PPN parameters, J 2 the Sun quadrupole moment. Re the Sun equatorial radius, a the semimajor 
axis of the orbit considered, e its eccentricity and n = {GM /the mean motion. In this analysis, we 
use 7=1 since this parameter is better determined by other types of observations like the Shapiro time 
delay (Bertotti ct al. 2003) and by Gaia itself that will be able to constrain it at the level of I0“® by observing 
light deflection (Mignard & Klioncr 2010). 

The fact that a large number of asteroids are considered with various different orbital parameters helps 
to decorrelate the two parameters. In this simulation, we obtain a sensitivity around aj^ ^ 10“^ and ap ~ 
7 X 10“^ with a correlation coefficient of 0.56 between the two parameters. These sensitivities are 1 order of 
magnitude lower than the ones obtained with planetary ephemerides (Pitjcva & Pitjev 2014; Vernia et al. 2014; 
Fienga ct al. 2015). The consideration of the complete set of asteroids (360,000 instead of 10,000) will improve 
slightly our current estimation. Planetary ephemerides will likely remain more powerful but Gaia will provide 
an independent and complementary estimation. 

In addition to these two parameters, we also consider a violation of the Strong Equivalence Principle (SEP). 
Such a violation appears in all alternative theories of gravitation. One effect produced by a violation of the 
SEP is that the trajectories of self-gravitating bodies depend on their gravitational self energy Pi (violation of 
the universality of free fall). It is characterized by a difference between the gravitational and the inertial mass 
usually parametrized by the Nordtvedt parameter p 

mg = mi + p^, (3.2) 

where mg is the gravitational mass and mi is the inertial mass. Using the same modelling as in (Mouret 2011), 
the estimated uncertainty on p using simulations of 10,000 asteroids is 9 x 10“"^. This parameter does not 
change the estimations on the fd and J 2 parameters and is not correlated to them. The only constraint available 
currently on p is provided by Lunar Laser Ranging measurements and is at the level of 4.5 x 10“^ (Williams ct al. 
2009). 

Moreover, in the PPN framework, the p parameter is unambiguously related to the PPN parameters (Will 
1993) by the relation 

,y = 4/3-7-3. (3.3) 

Instead of estimating 3 independent parameters (J 2 , fd and 7 ), one can introduce the last relation and estimate 
only two of them (J 2 and jd). By doing so, the estimated uncertainty on J 2 and jd becomes crjj ~ 9 x 10“® and 
( 7,3 ^ 2 X 10“^. Considering a violation of the SEP predicted by the PPN framework leads to an improvement in 
the estimation of the /3 PPN parameter. Moreover the correlation coefficient between (d and J 2 drops from 0.56 
to 0.18, which can be a considerable improvement. Similar conclusion holds for planetary ephemerides analysis 
and this gives a strong motivation to consider violation of the SEP with planetary ephemerides gravitation tests. 


3.2 Lense-Thirring effect 

The Lense-Thirring effect is a purely relativistic frame-dragging effect produced by the angular momentum of 
bodies. While the Lense-Thirring effect from the Earth has been detected with the LAGEOS satellite (Ciufolini & Pavlis 
2004) (see nevertheless the controversy raised in lorio ct al. (2011)) and with the Gravity Probe B mis¬ 
sion (Evcritt et al. 2011). The Sun Lense-Thirring has never been directly detected so far. Such a direct 
measurement would provide a new way to estimate the Sun angular momentum, an important quantity to 
assess interior models of the Sun (Pijpcrs 2006). Possibilities to measure the Sun Lense-Thirring effect has 
been mentioned in (lorio ct al. 2011; lorio 2011, 2012b, a). Nevertheless, as mentioned in (Folkncr ct al. 2014), 
planetary ephemerides analysis does not allow to disentangle the Sun Lense-Thirring from the Sun J 2 (even with 
the latest Messenger data). The two effects are completely correlated. This is partially due to the fact that all 
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planets are orbiting in a very similar plane with a nearly circular orbit. Considering the asteroids can help since 
they provide a wide range of different orbital parameters especially including larger inclinations that may reduce 
this correlation. With the 10,000 asteroids considered in this study, the uncertainty on the Sun Lense-Thirring 
is a factor 5 larger than its actual value. This means that a direct detection of this effect with Gaia observations 
seems difficult. This conclusion should be reassessed by the consideration of the full set of asteroids. Moreover, 
the combination of the Gaia dataset with radar observations performed at UGLA (Margot & Giorgini 2010) 
may also improve this conclusion. Finally, it would be interesting to assess the gain that an analysis combining 
planetary ephemerides with Gaia observations would bring. 

Nevertheless, even if a detection of the Sun Lense-Thirring seems to be unreachable, the fact to not include 
this effect in the modelling of the orbital dynamics produce biases in the estimation of the other parameters. 
Our simulations indicate that not including the Lense-Thirring in the equations of motion leads to a bias at the 
level of 10“® on J 2 and at the level of 5 x 10“^ on the (3 PPN parameter. Similar conclusions seem to hold for 
planetary ephemerides, which is also mentioned by lorio et al. (2011) and Folkner et al. (2014) for the J 2 . 

3.3 Fifth force formalism 

The fifth force formalism consists in a Yukawa modification of the Newtonian potential (Talmadge et al. 1988; 
Adclberger et al. 2009). This Yukawa modification is parametrized by a range of interaction A and by an 
intensity a. Gonstraints on these two parameters can be found for example in Fig. 31 of (Konopliv et al. 2011). 
At astronomical scales, the constraints are provided mainly by Lunar Laser Ranging and by ranging to Mars 
spacecraft. With simulation of 10,000 asteroids, we estimate the uncertainty on a for different fixed values of A. 
At the level of A = 10^° m, the uncertainty on a is slightly larger than the one from ranging to Mars spacecraft 
while for A = 10^^ m, our estimated uncertainty improves the current constraint by a factor 5. Therefore, this 
seems to be a very promising test. The correlation with the Sun mass needs nevertheless to be assessed (see 
the related discussion by Konopliv et al. (2011)). 

3.4 Temporal variation of the gravitational constant 

A lot of alternative theories of gravitation promotes the gravitation constant G to a dynamical field (typically to 
a scalar field). In this case, G becomes space-time dependant and can evolve for example with the cosmological 
evolution. Therefore, this class of models predicts a non-vanishing value of G in the Solar System. Our covariance 
study of the observations of 10,000 asteroids over 5 years indicates that Gaia will be able to constrain a linear 
evolution of the gravitational constant G/G at the level of 10“^^ yr~^ (or more precisely d\nGMQ/dt). The 
correlation with the Sun J 2 is very weak (see also (Mouret 2011)). The best current estimations on G/G are at 
the level of 10“^^ yr~^ and are provided by planetary ephemerides (Konopliv et al. 2011; Fienga et al. 2015). 

Moreover, very recently, it has been reported that the measurements of the gravitational constant G seem 
to undergo some periodic variations (Anderson et al. 2015) with a 5.9 years period. This analysis has been 
confirmed by Schlammingcr et al. (2015) even if the conclusion has been weakened. We have performed a 
simulation including a periodic variation of G. This shows that asteroid observations from Gaia will be able to 
constrain the relative amplitude of the oscillations at the level of 10“^°. This accuracy is 5 orders of magnitude 
smaller than the amplitude predicted by Anderson ct al. (2015); Schlamminger et al. (2015). No correlation is 
expected with the Sun J 2 . This means that Gaia will be able to rule out such a temporal variation of G. Note 
that, as mentioned by lorio (2015), the current planetary ephemerides analysis is also able to rule out this time 
variation with a similar accuracy. 

3.5 Standard Model Extension formalism 

The Standard Model Extension (SME) formalism has been developed to systematically describe violations of 
the Lorentz symmetry in all sectors of physics, including the gravitational sector. Gonsequences of SME on 
gravitational observations are developed into details in (Bailey & Kostelecky 2006; Kostelecky & Tasson 2011). 
In the linearized gravity approximation, the gravitational sector of the minimal SME is described by a space- 
time metric that depends on a symmetry trace-free tensor . It can be shown that the orbital dynamics is 
insensitive to (Bailey & Kostelecky 2006) and therefore depends on 8 independent parameters: , 

and . The heliocentric equations of motion can be found 

in (Bailey & Kostelecky 2006). 


Tests of GR with Gaia 


241 


The covariance analysis performed by considering 10,000 asteroids leads to the estimated uncertainties 
presented in Table 1. These uncertainties are better than the current best estimations of the SME parameters 
available in the literature (Kostelecky & Russell 2011). In particular, they are better than the ones obtained 
with planetary ephemerides (lorio 2012c; Hees et al. 2015). This is due to the variety of the asteroids orbital 
parameters while planetary ephemerides use only 8 planets with similar orbital parameters (same orbital planes 
and nearly circular orbits). Therefore, the estimation of the SME parameters with planetary ephemerides 
are degraded by these correlations (see the discussion in Hees ct al. (2015)). Using our set of asteroids, the 
correlation matrix for the SME parameters is very reasonable: the three most important correlation coefficients 
are 0.71, -0.68 and 0.46. All the other correlations are below 0.3. Moreover, the Sun J 2 is not correlated 
to these parameters. Therefore, Gaia offers a unique opportunity to constrain Lorentz violation through the 
SME formalism. A combined analysis with planetary ephemerides analysis (Hees et al. 2015), Lunar Laser 
Ranging (Battat et al. 2007) and atom interferometry (Miillcr ct al. 2008) would also be very interesting. Our 
analysis needs to be extended to include gravity-matter Lorentz violation parametrized by (d)'^ coefficients in 
the SME framework. 


Table 1. Sensitivity on the SME gravity parameters. 
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4 Conclusion 

In this communication, we have presented different possibilities to use asteroid observations from Gaia to 
perform different test of the gravitation theory. The estimation of the /3 PPN parameter and of a linear time 
dependance of the gravitational constant {d In GMq/ dt) would not be as good as the current estimations from 
planetary ephemerides. Nevertheless, the expected constraints are very complementary since they do not suffer 
from the same systematics errors. 

Moreover, we have shown that the asteroid orbits are sensitive to a violation of the Strong Equivalence 
Principle through the rj parameter. In the framework of the PPN formalism, this 77 parameter is related to the 
standard PPN parameters and this relation can help to reduce the correlation between the Sun J 2 and the PPN 
parameter (3. Similar conclusion holds for planetary ephemerides analysis. 

We have shown that the Sun Lense-Thirring effect is too weak to be detected with such dataset. Never¬ 
theless, it would be interesting to see if a combination of the Gaia dataset with radar observations of aster¬ 
oids (Margot & Giorgini 2010) or with planetary ephemerides can lead to a first dynamical detection of the 
Sun angular momentum through the Lense-Thirring. Nevertheless, we have shown that not including the Sun 
Lense-Thirring in the equations of motion leads to a bias in the estimation of the Sun J 2 and of the fd PPN 
parameter. Similar conclusions hold for planetary ephemerides analysis. 

The asteroid trajectories are also sensitive to a fifth force and can be used to improve constraints in this 
formalism. 

Finally, the number of asteroids and the variety of their orbital parameters provide a unique opportunity 
to constrain Lorentz violation through the Standard Model Extension formalism. The wide orbital parameters 
allows to decorrelate the SME parameters and will allow us to produce the best estimations on the SME 
parameters. 

Extended simulations considering the full set of asteroids (360,000 asteroids instead of 10,000 considered in 
this analysis) with refined astrometric precision (instead of 0.2 mas) and possible mission extension (6 years) are 
currently on-going. We are also currently assessing the improvement provided by combining the Gaia dataset 
with radar observations (Alargot & Giorgini 2010) that are complementary in the time frame and orthogonal 
to astrometric telescopic observations. 
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